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^ , Abstract 

The time evolution of driven two-level systems in the far off-resonance regime 
is studied analytically. We obtain a general first-order perturbative expression for 
^vq ■ the time-dependent density operator which is applicable regardless of the coupling 

| strength value. In the strong field regime, our perturbative expansion remains 

VO ■ valid even when the far off-resonance condition is not fulfilled. We find that, 

in the absence of dissipation, driven two-level systems exhibit coherent population 
trapping in a certain region of parameter space, a property which, in the particular 
case of a symmetric double-well potential, implies the well-known localization of 
the system in one of the two wells. Finally, we show how the high-order harmonic 
generation that this kind of systems display can be obtained as a straightforward 
■ application of our formulation. 



oo 



CO 

o 



a,: 
^— > 



PACS number(s): 42.50.-p, 42.50.Ct 



•rn . I. INTRODUCTION 

X 

Numerical investigations in driven symmetric bistable systems have revealed a num- 
ber of striking quantum phenomena Q. On one hand, Grossmann et al. || have shown 
that the parameters of the external driving field can be appropriately tuned so as to 
produce coherent suppression of tunneling, a property that can be used to localize the 
quantum system in one of the two wells. By properly choosing the shape of the driving 
laser pulse, localization can be achieved even if the system is initially in a delocalized 
eigenstate: the driving field can take it into a localized state and then keep it there ||. 
On the other hand, this system has also been shown to exhibit high-order harmonic 
generation in a region of parameter space which overlaps to a certain extent with that 
where localization occurs Q. Harmonic generation is a consequence of the fact that, 
in a strong field, the induced dipole moment responds with frequencies that are integer 
multiples of the laser frequency, thus giving rise to the appearance of the corresponding 
peaks in the emission spectrum. 
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An interesting aspect of driven double-well systems is that many of their relevant 
features can be captured in a simple two-level model. Using this kind of approach it 
has been shown || that, in a wide parameter range Q, the localization conditions can 
be correctly obtained as a zeroth-order result of a perturbative analysis in the small 
parameter Aq/wl, with Aq being the transition frequency of the two-level system and 
lul the driving field frequency. Specifically, in this far off-resonance regime (Ao/wl "C 
1), localization was found to occur at the zeros of the Bessel function Jo (2Qo/ w l)) 
where Slo is the Rabi frequency. In contrast, first order perturbation theory is required 
at least to account for the high-order harmonic generation that occurs in the strong 
field regime of driven two- level systems ]?]]. In this respect, several approaches have 
been developed, primarily aimed at obtaining a perturbative first-order solution for 
the equation of motion governing the evolution of the induced dipole moment |8|-1C]. 
All these approaches have to deal with secular terms appearing at first order. These 
terms become divergent in the long time limit and, consequently, have to be carefully 
resummed in order for the perturbative solution to be applicable at any time. 

A somewhat different approach, specifically designed for the strong field limit, has 
been followed in Ref. [11], where a first-order perturbative solution was obtained by 
using a series expansion dual to the Dyson series, and by resumming the corresponding 
secular terms by renormalization group methods [{12j|. 

Since localization can occur in a two-level system even when the strong field condi- 
tion is not satisfied, it is interesting to have a perturbative solution applicable regardless 
of the coupling strength value. Such a solution would permit to treat within a unified 
formulation the localization and harmonic generation properties exhibited by this kind 
of systems. On the other hand, and because of the extensive use of two-level systems 
as a first approximation to treat more complex physical systems, analytical solutions 
of this type become of particular interest. 

In this paper we derive a general first-order analytical solution for the perturbative 
evolution of a driven two-level system in the far off-resonance regime. Specifically, a 
first-order expression for the time-dependent density operator of the system, which is 
applicable regardless of the coupling strength value, is obtained. Remarkably, in the 
strong field limit, our perturbative expansion turns out to be valid even away from the 
far off-resonance condition, and therefore, in this particular case, it becomes indistin- 
guishable from that of Ref. [11|. Moreover, our zeroth-order Hamiltonian already in- 
cludes all the slowly varying contributions thus preventing, unlike previous approaches, 
secular terms from appearing in the corresponding first-order results. 

In Sec. II we develop the formalism and obtain the dynamical evolution of the 
corresponding density operator. Then, in Sees. Ill and IV this formulation is applied, 
respectively, to the study of localization and high-order harmonic generation. We find 
that, in the absence of dissipation, driven two-level systems exhibit coherent population 
trapping in the far off-resonance regime. In the particular case of a symmetric double- 
well potential such a property implies the well-known localization of the system in one 
of the two wells. Finally, the main conclusions are summarized in Sec. V. 
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II. TIME-DEPENDENT DENSITY OPERATOR 



We consider a two-level system driven by a linearly polarized laser field of frequency 
cjl and amplitude Eo- The energy difference between the upper level state |2) and the 
lower level state |1) is denoted Ao- In the dipole approximation the Hamiltonian reads 
(in atomic units, h = 1) 

Hd = ~Y (0"22 - CTil) - VLq 9 (t) COS (u L t) ((712 + 0"2l) , (1) 

where <7jj = \i){j\ is the transition operator and Qo = Eofi is the Rabi frequency, with 
\x being the (real) dipole matrix element between |1) and |2). It has been assumed 
that the laser polarization vector and the dipole moment are oriented along the same 
direction. Moreover, we have included in Eq. (|l|) an envelope function g (i) in order to 
account for slow variations of the intensity, as occurs at the turn on and off of the laser 
pulse. 

Appropriate scaling of the time-dependent Schrodinger equation permits the identi- 
fication of the relevant dimensionless parameters. In this respect, especially suitable for 
our purposes is the variable change r = uj-^t, which yields the following dimensionless 
Hamiltonian 

H = —5- (022 - o-n ) - — / (r) cos (r) (an + <t 2 i) , (2) 
2w L 

where / (r) = g(r/uj L ). 

On the other hand, particularly convenient for studying the localization properties 
associated with the above Hamiltonian are the following coherent superposition states 

| r ) = -L(|l) + |2», (3) 

|1> = ^(|1>-|2», (4) 

which, in the case of a symmetric double-well potential, correspond to states localized 
on the right and left wells, respectively. In the basis {|r), |1)} the Hamiltonian (^) clearly 
displays the symmetry of the system. Indeed, in this basis Eq. (||) takes the form 

H = -—9- (air + o- r i) - — / (r) cos (t) (cr rr - an) , (5) 
2lu l lul 

which, for a driving field with constant amplitude (/ (r) = 1), is manifestly symmetric 

under the combined transformation r < > 1 and r — > r + tt. 

As already said in the Introduction, in the far off-resonance regime, i.e., when 
Ao/wl <C 1, localization occurs (for / (r) = 1) at the zeros of the Bessel function 
Jo (2£Io/ljl). This is the regime in which we will be primarily interested. Thus, Ao/cjl 
is going to be, in principle, the small parameter of our perturbative analysis, which, 
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consequently, will be applicable regardless of the value of the coupling strength On/wL- 
That is, in the far off-resonance regime our perturbative series will remain valid in both 
the strong and the weak field limit. 

In order to obtain a small Hamiltonian, proportional to Ao/wl, we perform the 
following unitary transformation 



U(r) 



-i0(r)(cr rr -CT 11 ) 



with 



/dr—f (r) cos (r) ~ — / (r) sin (r) 



(6) 



(7) 



where, in the last step, we have used the fact that the pulse envelope / (r) is a very 
slowly varying function within a laser period. The transformed Hamiltonian now takes 
the form 



H' = UHU^ 



iUU A 



A, 



2^ L 



(8) 



where U + denotes the derivative of the adjoint of U with respect to the dimensionless 
time r. 

In order to prevent secular terms from appearing at first order in the perturbative 
expansion, it is necessary to incorporate all the slowly varying terms in the zeroth-order 
Hamiltonian. With this purpose, and using again the fact that / (r) hardly changes 
in a laser period, we express the time-dependent coefficients in Eq. (||) as the Fourier 
series 



where 



3 ±2i0(r) 



A + A±(r), 



An 



+oo 



Jo 

'2n 



2^o 



f(r) 



f(r) 



A±(r) = £(±1)V» 

n=l 

Substituting Eq. (||) into Eq. (g) we then arrive at 

H' = Hq + AH'. 



e inr + (_i)» e - 



with 



An_ 

2^ L 



Ao (o"ir + er ri ) 



(9) 

(10) 
(11) 

(12) 
(13) 
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AH' = -p- (A + (r)a h + A_(r)a rl ) . (14) 

The slowly varying part H' Q is going to be considered our zeroth-order Hamiltonian 
while, in the far off-resonance regime, AH' becomes a small perturbation. Incidentally, 
note that in the strong field limit, that is, when 

C = — -/ (t) ^> 1, for any r (15) 



the Bessel functions entering Eq. ( |Tl| ) become, for n < n c ~ ( [1 



while, for n > n c they decay very fast as J n (C) ~ (eQ/2n) n /\/2irn. Hence, the pertur- 
bation AH' takes now the form 

AH' = — A ° (U + (r)a h + n_(r)a rl ) , (17) 

2^uj L n f{T) 

with 

n±(r) » (1/V5F) £>1)» (e^ + (-l)"e— ) cos (C - f - • (18) 

n=l 

Thus, in the strong field limit, AH' turns out to be proportional to the parameter 
Ao / y^L^o/ ( r )> which has the interesting consequence that whenever the latter be- 
comes small (for any r) our perturbative results remain valid irrespective of the value 
of A /uj l . 

We now proceed to solve perturbatively the quantum Liouville equation for the 
transformed density operator p' , which contains all the dynamical information about 
the system, 

% = -i [H', p'} = C!p' = (C' + AC) p>, (19) 
P '(t) = U(t)p(t)U + (t). (20) 



The linear operators C and AC' entering Eq. (|19|) represent the quantum Liouville 
operators corresponding to the Hamiltonians Hq and AH' , respectively. This evolution 
equation can be exactly solved to the lowest order. In doing so, one finds the following 
expressions for the matrix elements of the density operator in the states |1) and |2), 

/4 0) '(r)=/4(0), i = l,2 (21) 
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p ( S'(r)=p' 12 (0)e~ T . (22) 



It is worth noting that according to Eqs. (|20| ) and fl2l|), to the lowest order, the 
populations of the states U + (t)\1) and U + (r)\2) remain constant. This is a direct 
consequence of the fact that, for / (r) cte, the states 

U+(t)\1) =cos0(r)|l) +ifrin0(r)|2), (23) 

U + {t)\2) =sin^(r)H) - icos0(r)|2), (24) 

become the zeroth-order Floquet states of the system corresponding to the quasienergies 
— AoAo/2 and +AoAo/2, respectively. 

The zeroth-order density operator p(°)' can be written as 

,(°)'(r) = e^V(0) = E + (p / i 2 (0)e^ T |l)(2| +h.c.) . (25) 

i=l,2 ^ ' 

Notice that this expression provides, in turn, an useful operational definition for the 
evolution operator e^o T . 

With the aim of deriving a first-order perturbative expression for the time-dependent 
density operator, we rewrite the quantum Liouville equation fll9|) in the interaction rep- 
resentation with respect to C' . Specifically, 

^=AC[p' l (r), (26) 

where 

pi = e- c ^p\r), (27) 

AC[ = e - c 'o T A£'e c '° T . (28) 
The solution of Eq. (pq ) can be expressed as the following infinite series 

p' l (r) = p'M+ f dnA£[(T 1 )p[(0)+ f dn r dT 2 AC[(T 1 )AC[(r 2 )p[(0) + ... (29) 
Jo Jo Jo 

This is a perturbative expansion in the small parameter Aq/oj^ which, in the strong 
field limit, becomes an expansion in the small parameter Ao/y^L^o/ (r), as already 
said. Thus, up to first order, one obtains 

p'{r) = e £ oy (0) + e c '" T f dne^ Tl A£'(n)e £ ° V(0) + O (e 2 ) , (30) 

Jo 

where e represents, in general, the small parameter characterizing the expansion. On 
the other hand, the integrand on the right hand side of Eq. (|30|) takes the form 
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= E &fa)|i>(<l + (ei 2 (ri)e i ^r I(r - ri) |l)(2| +h.c. 

i=l,2 ^ 



(31) 



where 



e' fe (r) = -z(i| Aff'(r),e £ i>'(0) 



(32) 



and the definition of the propagator e c '^ T Tl \ as provided by Eq. (p5|), has been used. 
Substituting then Eq. (p5|) into Eq. (p2]) one obtains, after some lengthy algebra, 



+oo 



4-(r) = (-1) J '^E j2 «+! (0 sin t( 2n + ^ ( / 5 i2(0)e'~ T + h.c. ) , (33) 



n=0 



A 



&(t) = E J 2«+i (0 sin [(2n + l)r] [p 22 (0) - p' n (0)] 



ra=0 



An ^ , ; A A T 

+ 2i_° V J 2 „ C cos (2nr p' 12 e «l r , 
^^i 



(34) 

£2i(t)=&(t), (35) 
with C = 2fin/ (r) A^l- Substitution of the above expressions into Eq. ( |3i"| ) and then 
into Eq. (pOI) leads to the following final result 



At) 



+ 



p'i 2 (0) + — /3(r) 



P22(0) 



Ac 



-air 



|2)(2| 



e^|l)(2| + h.c.^+0( e 2 ), 



(36) 



where 



a(r) ^ig-wi(c) f , 2(0) ^H + ^] T+e -*^)-^]^_ 2 ] +hx .| 

2 n =o + 1 [ \ / J 



(37) 



« - iE%#[^(°)-^(o)]^ 

n=0 ^ 



+ 2.E 



n=l 



^p' 12 (0)sin(2nr) 



i[(2n+l)-^]r + e -i[ { 2n+l)+^]r _ ^ 

(38) 
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Thus, the first-order density operator of the driven two-level system is finally given by 

p(r) = U + (t) P '(t)U(t), (39) 

with U(t) defined in Eqs. (|6|) and (|7|). This expression, which contains no secular 
terms, is our central result and is applicable not only in the far off-resonance limit, i.e., 

AoM « 1, (40) 

but also in the regime where 

Aq/wl « V^of (r) /w L » 1. (41) 



III. COHERENT POPULATION TRAPPING AND 
LOCALIZATION 

Coherent population trapping and localization properties are conveniently analyzed 
in terms of the diagonal matrix elements of p(r) in the basis {|r), |1)}. By using Eqs. 
(H) and (|39|) one can express the population (r|p(r)|r) as 

pAt)= 1 -{1+[p' 12 (t)+p' 21 (t)]}. (42) 

Taking into account that p'(0) = p(0) one obtains, after substitution of Eq. ( p2|) into 
Eq. (fi"2^), the following zeroth-order expression 

Pi°r ) (r) = l + \pu(0)\ cos (^°T + tp y (43) 

where pi2(0) = |/Oi2(0)|e l</3 . Equation (]43|) shows that, in the absence of dissipation, 
the populations of the |r) and |1) states, which in general oscillate with a frequency 
AoAo (curves (a) and (b) in Fig. 1), become trapped whenever the parameters of the 
driving field are so tuned that the Bessel function Aq = Jo (20q/ (t) /o>l) attains a zero 
(curve (c) in Fig. 1). Remarkably, such a population trapping occurs regardless of the 
system initial state. This is a general zeroth-order result that remains valid as long as 
conditions ( f40| ) or ( |4l| ) are satisfied. In the particular case of a symmetric double-well 
potential, |r) and |1) become states localized on the right and left wells, respectively, so 
that, when the system is initially in either of the two wells, one obtains the peculiar 
localization property reported in previous works f|-||]. 

Substituting Eq. fl36| ) into Eq. (|42| ) one finds the following final expression for the 
first-order time-dependent population 
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Figure 1: Zeroth-order population p T i(r) for a pulse envelope /(r) = 1 — exp(— r/10) and 
Ao/cjl = 0.2, |p 12 (0)| = y/3/4, <p = tt/4; (a) ft„/wL = 0; (b) Q Q /to L = 1.8; (c) Sl /w L = 1.202. 



/M r ) 



- + Pi2(0) cos r + ip 

2 \ 



+ 



[p 22 (0)-pn(0)]^ 



n=0 



■W (C) ( ,„ , u AqAq 

cos( In + l)r — cos r 

2n + l I 



2^|p 12 ( 



y^) sin (2nr)sinf^r + ^+0( e 2 ). (44) 
i 2n \ / 

n=l ' 



Equation (44) shows that now the populations evolve undergoing rapidly oscillat- 
ing changes of the order of e, superimposed to the much slower dominant oscillating 
evolution of frequency AoAo (see curve (a) in Fig. 2). In general, these high frequency 
oscillations cannot be eliminated by the external field and as a consequence, unlike the 
previous case, to this order it is not possible to achieve exact coherent trapping for 
any system initial state. Indeed, when A = the population p vv {r) evolves describing 
small-amplitude rapid oscillations (with frequencies that are integer multiples of the 
laser frequency) about the initial population p rr (0) [curve (b) in Fig. 2]. Yet, exact 
first-order coherent trapping can still be obtained by properly choosing the initial state. 
Taking pn(0) = p 2 2(0) = 1/2 and ip = 0, Eq. (H) becomes, when A = 0, 



Prr(T) = - + pn(0) = Prr(O). 



(45) 
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Figure 2: First-order population /3 rr (r) for /(r) = 1 and A /wl = 0.2; (a) pn(0) = 3/4, 
p 22 (0) = 1/4, |p 12 (0)| = V3/4, ^ = tt/4, Oo/^l = 1.8; (b) p u (0) = 3/4, p 22 (0) = 1/4, 
|pia(0)| = V3/4, y> = tt/4, n /^ L - 1.202; (c) p n (0) = p 22 (0) = 1/2, |p 12 (0)| = 1/2, ip = 0, 
^o/cjl = 1.202. 

That is, when the system is initially in either the state |r) or |1), first-order coherent 
population trapping can be achieved by simply choosing the parameters of the driving 
field in such a way that Ao = (curve (c) in Fig. 2). Notice that, as before, in the 
special case of a symmetric double-well potential, this result implies the localization of 
the system in one of the two wells. 

On the other hand, it is interesting to note that according to Eq. ([4iD the population 
of a system initially located in the right well of a symmetric double-well potential, 
evolves in the limit of zero driving field as 

Prr (t) = cos 2 (^-tj (46) 

so that, the well-known result for the undriven system is recovered, as it should be. 

Notice that the pulse envelope f(t) simply modifies the way in which the system 
reaches its steady state (see Fig. 1). As a consequence, localization cannot be pertur- 
batively achieved starting from an initial delocalized eigenstate. Indeed, as is apparent 
from Eq. (|44|), in this case the system remains always delocalized, within terms of the 
order of e. This result, which is in contradiction with the numerical results obtained in 
Refs. HQ, simply reflects that such a localization process is a nonperturbative effect 
that cannot be accounted for by our treatment. 
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IV. HARMONIC GENERATION 



Driven two-level systems has been extensively used as convenient models to under- 
stand the basic mechanism underlying the phenomenon of high-order harmonic gen- 
eration IPI^.IC, 14-1^1- This is so primarily because, as pointed out by Sundaram and 
Milonni ]7|, such a simple system already exhibits the main features observed exper- 
imentally in the emission spectrum of atoms in very intense laser fields, namely, the 
existence of a plateau in the harmonic spectrum followed by a sharp cutoff. This, in 
turn, is an indication of the fact that such features are intrinsic properties of strongly 
driven systems. 

The purpose of this Section is to show how the previous formulation leads straight- 
forwardly to the well-known results of high-order harmonic generation in strongly driven 
two- level systems. 

The coherent part of the emission spectrum S(uj) is proportional to |(i(u;)| 2 , with 
d(cj) being the Fourier component at the frequency u of the induced dipole moment 

(dm an 



to+T 



dte tuJt (d(t)) 



to 



(47) 



The time evolution of {d(t)) follows directly from Eq. 
that the quantum dipole operator is given by 



f). Indeed, taking into account 



d = p(a 12 +<7 2 i) = At (cr r 



(48) 



one finds 



(d(i)> = Tr [p{t)d] = p [p„{t) - p n (t)} = 2p (p TT (t) - 1/2) , (49) 

where, in the last step, we have used that p rr (t) + p\\(t) = 1. Therefore, by substituting 
Eq. ( f44|) into Eq. (fi"9|), and then this latter into Eq. ( [47|) one can immediately obtain 
the corresponding emission spectrum. In doing so, one finds the well-known result that 
S(u) consists, in general, of three types of peaks: 

i) Low frequency components located at u = ±AoAo and intensities proportional 
to 

,2_..2,„ (50) 



l rf ±A Ao| 



ii) Hyper-Raman lines located at frequencies 2nui\ J ± AoAo and intensities propor- 
tional to 



\d 



2najL±A()Ao I 



/i P12(0) 



2n 



n 



1,2,3,... 



(51) 



iii) Odd harmonic components located at frequencies (2n + 1) oj^ and intensities pro- 
portional to 
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d 



(2n+l)oj L | 



(J>— [P22(0) 



Pn(0)] 



J2n+1 (C) 



2n+ 1 



7? 



0,1,2, 



(52) 



Formulas of this type have been previously obtained by Ivanov and Corkum || and 



Dakhnovskii and Bavli [1C|. 

It is interesting to note that exact first-order coherent population trapping (local- 
ization) has a twofold manifestation in the system spectroscopic properties: a zero 
frequency peak appears and the spectrum consists of only even harmonic components. 

On the other hand, the appearance in the strong field regime (£ 3> 1) of the plateau 
and corresponding cutoff in the harmonic spectrum follows from the asymptotic be- 
haviour of the Bessel functions, which for n < n c ~ £ are given by Eq. ([!(]) while for 
n > 



n, 



they drop very rapidly as (e(/2n) n /^/2~Kn. 



V. CONCLUSION 

Driven two-level systems display a number of interesting features that are common 
to more general driven systems. This fact, makes them convenient starting points for 
analyzing more complex physical systems. Moreover, under certain circumstances, they 
allow for an analytical treatment which would be otherwise impossible. Such analytical 
treatments provide detailed information on the system response to specific changes in 
the external parameters, and this information can be used to control its dynamical 
evolution. 

In this work we have studied analytically the time evolution of driven two-level 
systems in the far off-resonance regime (Ao/^l *C 1). This was done by performing a 
perturbative analysis based on a convenient zeroth-order Hamiltonian which takes care 
of divergent secular terms. In this way, we obtained a general first-order expression for 
the time-dependent density operator which is valid regardless of the coupling strength 
value. Interestingly enough, in the strong field regime, this expression turns out to 
be applicable even when the far off-resonance condition is not satisfied. Indeed, the 
perturbative analysis remains valid as long as the condition Ao /cut <C y/Qof (t) /cjl 3> 
1 holds. 

The analytical formulation presented in this paper makes it possible to treat in 
a unified framework different aspects of the dynamical evolution of driven two-level 
systems. In particular, from the time evolution of the populations it immediately follows 
that the well-known phenomenon of tunneling suppression in a symmetric double-well 
potential can be considered as a specific manifestation of a more general population 
trapping phenomenon. 

To the lowest order, regardless of the system initial state, the populations of the 
coherent superpositions |r) and |1) become trapped, in the absence of dissipation, when- 
ever the parameters of the driving field are tuned in such a way that the Bessel function 
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Jo (2fio/ (t) /wl) vanishes. To first-order, however, small-amplitude rapid oscillations 
in the time evolution of the populations appear, which, in general, cannot be eliminated 
by the external field. As a consequence, to this order, such a result only remains valid 
in an approximate way. Nonetheless, exact first-order coherent population trapping 
still occurs when the system is initially in either the state |r) or |1). In the particular 
case of a symmetric double- well potential, in which |r) and |1) become states localized 
on the right and left wells, respectively, this result implies the localization of the system 
in one of the two wells. 

The present formulation also leads straightforwardly to the well-known results of 
high-order harmonic generation in strongly driven two-level systems. 

Of course physical systems exhibit a number of interesting features which are of 
nonperturbative origin and, consequently, cannot be accounted for by a perturbative 
treatment of the type presented here. In spite of this, the present formulation can 
be useful as a starting point for analyzing a large variety of driven systems in the 
perturbative regime. 
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